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r~H | Supernovae (SNe) driven winds are widely thought to be very influential in the 

O t high-redshift Universe, shaping the properties of the circum-galactic medium, enrich- 

ing the intergalactic medium (IGM) with metals and driving the evolution of low-mass 
galaxies. However, it is not yet fully understood how SNe driven winds interact with 
their surroundings in a cosmological context, nor is it clear whether they are able to 
significantly impact the evolution of low-mass galaxies from which they originate by 
altering the amount of cold material these accrete from the cosmic web. Indeed, due 
to the strong constraints on resolution imposed by limited computational power, all 
cosmological hydrodynamics simulations to date resort to implementing more or less 
physically well motivated and complex subgrid models to trigger galactic winds. To 
explore this issue, we implement a standard Taylor-Sedov type solution, widely used 
in the community to depict the combined action of many SN explosions, in a cosmo- 
logical resimulation of a low mass galaxy at z ^ 9 from the 'nut' suite. However, in 
contrast with previous work, we achieve a resolution high enough to capture individual 
SN remnants in the Taylor-Sedov phase, for which the Taylor-Sedov solution actually 
provides an accurate description of the expansion. We report the development of a 
high- velocity, far-reaching galactic wind produced by the combined action of SNe in 
the main galaxy and its satellites, which are located in the same or a neighbouring 
dark matter halo. Despite this, we find that (i) this wind carries out very little mass 
(the measured outflow is of the order of a tenth of the inflow/star formation rate) 
and (ii) the cold gas inflow rate remains essentially unchanged from the run without 
SNe feedback. Moreover, there are epochs during which star formation is enhanced 
in the feedback run relative to its radiative cooling only counterpart. We attribute 
this 'positive' feedback to the metal enrichment that is present only in the former. 
We conclude that at very high redshift, efficient SNe feedback can drive large-scale 
galactic winds but does not prevent massive cold gas inflow from fuelling galaxies, 
resulting in long-lived episodes of intense star formation. 

Key words: methodsmumerical-galaxies: high-redshift-galaxies: formation-galaxies: 
evolution-intergalactic medium-super novae: general 



1 INTRODUCTION 

Significant outflows of gas from star-forming galaxies, known 
as galac tic winds, have been observed both in lo cal galax- 
ies (e.g. iLehnert fe Heckmanl [1995 iMartinl [l999h and high 

redshift (z ~ 3) Lyman break gala xies (LBGs) (e.g. 

IPettini et all l200ll : IShaplev et all 1 20031 ). These winds are 
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usually associated with starbursting galaxies (for example 
LBGs have star formation rates (SFRs) of ~ 100M©/yr) 
and in this case are often referred to as 'superwinds'. These 
winds can reach large distances from the source galaxy and 
are typically observed to have temperatures greater than the 
'escape temperature' of the potential well of t heir host halo 
e.g. M82, a classic example of a superwind (iLehnert et al.l 
1999). This suggests that at least some of their metal en- 
riched gas may be able to escape into the IGM. 
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The production of these superwinds is attributed to si- 
multaneous SNe explosions whose remnants can overlap pro- 
ducing a bubble of hot, low densi ty gas that can 'blowout ' 
of the galaxy and become a wind (jMcKee fe Ostrikerl[l977h . 
The main requirement for 'blowout' is that the rate of 
SNe is high enough, such that the remnants overlap before 
they can cool radiatively (|Heckman et al.lll990l : I David et al.l 
Il990h . Observations of outflows from M82 are consistent 
with multiple SNe s upplying energy to drive the wind 
(|Heckman et al.|[l990h . Alternatively, it has been advocated 
that multiple coherent SNe, as seen in the superwind sce- 
nario, are not the only mechanism able to produce outflows 
from dwarf galaxies. Dwarf galaxies can lose a significant 
fraction of their mass on timescales of a Gyr, in a more qui- 
escent manner. Fuelled only by an average SFR, cold gas 
clouds in the interstellar medium (ISM) are evaporated by 
SNe remnants and form a w ind if the evapora ted gas reaches 
high enough temperatures (jEfstathioul l2Q00h . 

Due to their ability to eject material from galaxies, 
galactic winds are considered a potential solution to sev- 
eral outstanding problems in the field of structure forma- 
tion. They are proposed to play a role in the pollution 
of the IGM and intracluster medium with metals, the oc- 
curre nce of metallicity gra dients within individual galaxies 
(e.g. iHeckman et al.1 l2000h . the properties of dwarf galax- 
ies (Dek el fe Silkl 1 19861 ) and even to provide a reservoir 
of hot gas for re- accretio n by galaxies at lower redshifts 
(jOppenheimer et al-llioioh . 

The most influential role of galacti c winds, however, 
i s arg uably their impact on the IGM. iFurlanetto &; Loebl 
(2003) demonstrate analytically that LBGs can propa- 
gate ~ lOOkpc into the IGM, which is compatible with 

the sizes of observed Hi-deficient regions around LBGs 
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(Adelberger et al. 2003). Based on cosmological simulations, 
lAguirre et al.l (|200ll ) argue that SNe driven winds from mas- 
sive galaxies alone (Mbaryon > 1O 1O-5 M0) could pollute 
the entire IGM to its observed level. It has also been pro- 
posed, however, that protogalaxies (with masses comparable 
to present day dwarf galaxies) could be better candidates 
for polluting the IGM, since the lower velocity winds they 
would produce at z = 9 would not per turb the IGM in th e 
way that superwinds at z = 3 would (|Madau et al.ll200ll ). 
There is also evidence for metals in the IGM by 2; = 5, 
requiring pollution at an earlier epoch than that probed 
by observations of LBGs. S imulations of an isolated dwarf 
bv iMac Low &; Ferraral (1999) show that almost all metals 
are ejected from dwarf galaxies (with M gas < 1O 9 M0), even 
though the ejection of gas mass is surprisingly inefficient. 

It is clear that galactic winds have the potential to im- 
pact vast regions of space and influence structure formation 
on large scales. Therefore, the importance of studying this 
phenomenon in its full cosmological context cannot be over- 
estimated. However, while there have been many success- 
ful numerical si mulation studies of galac t ic winds in indi- 
vidual galaxies (IMac Low fe Ferraral 1 19991 : iMori et al.ll2002l : 
IScannapieco &; Briiggen 120101 ). capturing the production of 
a galactic wind in cosmological simulations remains a chal- 
lenge. In such simulations, galactic winds tend to be im- 
posed directly, rather than allowed t o develop naturally as 
a result of overlapping SNe bubbles (|Springel &; Hernauistl 
120031 : IScannapieco et al.ll2006l : iQppenheimer fe Davell2008h . 
We recall that an important ingredient of analytical stud- 



ies of SNe driven winds is a multiphase ISM i.e. a hot, 
rarefied atmosphere, containing cold den se gas clouds (e.g. 
iMcKee Ostrikerlll977l : lEfetathioul koOO). If the resolution 
of a cosmological simulation is not sufficient, the multiphase 
nature of the ISM is not captured and it is not possible to ac- 
curately model SNe remnants. In fact, the inclusion of a mul- 
tiphase ISM ' by hand' is a common feature of galactic wind 
mode ls (e.g. ISpringel &; Hernquistll2003l : IScannapieco et aD 
l2QQ6h . 

Examining galactic winds in a cosmological context, 
rather than in isolation, introduces further complications. 
Galaxies are now embedded in the cosmic web and are sub- 
ject to inflows as well as outflows. The mode of inflow will be 
dependent on the galaxy halo mass; low-mass galaxy haloes 
(Mvir < M s hock ~ 1O 12 M ) cannot host a stable shock at 
the virial radius and so primarily accrete cold (T < T v ir ) gaS 
whereas high- mass galaxy haloes (M v i r > M s h oc k) primarily 
accre te shock- heated, hot (T ~ T v i r ) gas (|Birnboim &; Dekel 
2003). It has been suggested that this framework, when cou- 
pled with active galactic nuclei (AGN) and SNe feedback, 
could naturally give rise to some of the bimodality in galaxy 
properties we observe but usually fail to reproduce with 
simulations (|Dekel &; Birnboiml l2006). Indeed, semi-analytic 
modelling has shown that incorporating such processes re- 
sults in an improved fit to most galaxy prope rty trends e.g. 
colou r bimodality, l umin osity function etc ([Cattaneo et al.l 
2006). iKeres et aD ([20091 ) have shown, however, that AGN 
radio-mode feedback (which can prevent hot gas from being 
accreted) has little effect on high-mass galaxies in cosmolog- 
ical simulations because most of these were built hierarchi- 
cally from lower mass galaxies which gained their baryons 
via cold accretion. They suggest it is feedback in low mass 
galaxies, such as SNe driven winds, that is of greatest impor- 
tance. It seems that studying accretion and feedback in low 
mass protogalaxies is, therefore, vital for our understanding 
of galaxy evolution as a whole. 

Recent cosmological simulations are in agreement that 
for low mass galaxie s (with M v i r < M s h oc k) most inflow- 
ing gas is cold (e g. IKeres et all 120051 : lOcvirk et aD 120081 : 
iBrooks et al.l l2QQ9h . All of these studies also demonstrate 
qualitatively that much of this cold inflow o ccurs along 
filaments (e g Fig 17 of IKeres et al.1 (l2005h. F i g. 5 of 
lOcvirk et all (|2008h and Fig. 5 of IBrooks et al.1 (|2009h ). 
What has not been quantified, however, is how much cold 
accretion occurs in a spherically symmetric way i.e. with 
the same geometry as assumed for hot accretion. While it is 
reasonable to assume that this spherical cold accretion can 
be safely ignored when the galaxy is surrounded by a halo 
of shock-heated gas (i.e. it has M v i r ~ M s h oc k), this is not 
the case when the halo is below the shock mass threshold 
and there is not necessarily anything to impede the spher- 
ical accretion of cold material from the IGM. Indeed, most 
of the visual demonstrations that cold accretion takes place 
along filaments in these previous studies focus on examples 
when the galaxy has a hot gas halo and cold accretion is nec- 
essarily confined to the filaments anyway. This paper thus 
complements such studies by measuring how cold gas accre- 
tion is divided between the filamentary and diffuse (spheri- 
cally symmetric) components in low mass (M v i r < M s h oc k) 
haloes, in simulations both with and without SNe feedback. 

Most of the aforementioned studies also incorporate 
SNe feedback, although it is unclear whether a substantial 
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Figure 1. Maps showing a z-projection of the hierarchy of grids (first row) and corresponding gas density (second row) in the feedback 
run at z = 9. From left to right: The whole box, a cube of side 12r v i r « 60kpc centred on the main halo and a cube of side 2r v i r also 
centred on the main halo respectively. The white circles in the middle and right gas density images (second row) indicate the virial radius 
of the main halo and its subhaloes. In the top left panel showing the whole box, we can see the 128 3 grid covering the whole image and 
the 3 nested grids of equivalent resolution 256 3 (orange), 512 3 (yellow) and 1024 3 (green) centred on our main halo. Additional grids 
can also be seen (blue, pink) indicating where additional levels of AMR have been triggered in the 1024 3 grid. The middle and right 
panels show the triggering of additional levels of AMR (up to a maximum of 15) in dense structures, such as the main galaxy disc and 
subhaloes. 



galactic wind ever develops and they do not address the spe- 
cifi c impact of this on t he ac cretion processes. Recent work 
by Ivan de Voort et al.l (|2010h attempts to tackle this issue 
with a suite of cosmological simulations employing different 
physics. We note, however, that their SNe feedback is im- 
plemented by giving gas particles a velocity kick and so the 
galactic wind is an input rather than a result that arises nat- 
urally. They conclude that the relative importance of cold 
accretion onto the galaxy halo is robust to changes to the 
feedback recipes, but that the relative importance of cold 
accretion onto the galaxy itself is decreased in the presence 
of SNe-driven winds. We revisit this question in this paper, 
but at much higher resolution and redshift (z ^ 9), where 
SNe blastwaves are individually resolved and so it is not nec- 
essary to put in a galactic wind by hand. We also emphasize 
that our high spatial and mass resolution allow us to fully 
resolve the filaments at very high redshift (z ^ 9), contrary 
to previous studies in which the spatial resolution, in partic- 
ular, is lacking in this redshift range. These resolution effects 
potentially enhance the likelihood of the filaments being de- 
stroyed by the SNe feedback, changing the balance between 
inflow and outflow and thus the subsequent evolution of the 
galaxy. 



A few high resolution studies of cosmological accre- 
tion in individual galaxies have been performed, but these 
tend to focus either on more massive g alaxies (i.e. with 
M vir > Mshock) at lo wer redshift (z ~ 3) (lAgertz et alf 2009: 
ICeverino et al.ll2010h or on the early formation stages of the 
first galaxies (M v i r ~ 5 x 1O 7 M ) at z > 10 where direct 
accretion of gas from the IGM b riefly precedes cold fila- 
mentary inflows (|Greif et al.ll2008h . We are primarily inter- 
ested in studying the epoch bracketed by these studies; that 
in which filamentary accretion is important for low mass 
haloes, during which galactic winds may also develop. 

In summary, for a low mass protogalaxy (i.e. with 
Mvir < M sri ock) at high redshift (~ z = 9 — 10), there are 
two main predictions about its interaction with the cosmic 
web : 1) a galactic wind will develop and extend far into 
the IGM, polluting it with metals and 2) cold gas will flow 
rapidly into the host halo via the web's filaments. What 
happens when these inflow and outflow mechanisms occur si- 
multaneously? We tackle this question with a suite of cosmo- 
logical resimulations (incorporating different physics) with 
0.5pc physical resolution in the densest regions, of a galaxy 
with Mvir ~ 5 x 1O 9 M at z = 9. At such high resolution, 
we can model individual SNe with a Sedov blastwave al- 
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lowing a galactic wind to arise naturally. We examine the 
mechanism via which the wind develops and measure its 
properties. Furthermore, by measuring the inflow and out- 
flow rates and comparing these with a control run (without 
SNe) we investigate whether the presence of a hot galactic 
wind can alter the accretion processes in a protogalaxy in 
such a way as to significantly impact its evolution. 

This paper is structured as follows. In Section [2] we de- 
scribe the suite of simulations we have performed and outline 
the gas physics employed (cooling, star formation, SNe feed- 
back) , further details of which are provided in the Appendix 
( Sect ion |5J). In Section[3]we distinguish between different gas 
phases in both runs, relating these to physical structures e.g. 
the filaments, satellite galaxies, the hot galactic wind (feed- 
back run only). In Section 2] we explore the mechanism re- 
sponsible for the development of the galactic wind in the 
feedback run and measure the wind properties. We measure 
the inflow and outflow in all gas phases in Section [5] and 
then relate the dominant accretion modes and outflows to 
the star formation rate in Section [6] comparing the cooling 
and feedback runs in order to isolate the impact of the SNe 
feedback on the formation and evolution of this galaxy at 
high-redshift. Finally, we outline our conclusions in Section 

m 



2 THE NUT SIMULATIONS 

The nut Simulations Q are a suite of ultra-high resolu- 
tion (maximum physical spatial resolution of ~ 0.5pc at all 
times) cosmological resimulations of a Milky Way like galaxy 
forming at the intersection of 3 filaments in a ACDM cosmol- 
ogy. To date we have performed 3 simulati ons, with differ ent 
gas physics, with the AMR code, RAMSES (Tevssier 200 j): i) 
adiabatic with a uniform UV background (|Haardt Madaul 
1996) turned on instantaneously at z=8.5 ii) cooling, star 
formation and UV background and iii) as ii) but with the 
addition of supernova (SN) feedback, including metal en- 
richment and its subsequent changes to the cooling func- 
tion. Simulations i) and ii) are presented elsewhere (Powell 
et al, in prep); for the sake of completeness we summarise 
the technical details below. 

The galaxy resides in a halo with a virial mass of 
1.05xlO n M© at z ~ 3 and has been selected such that it 
lies in an isolated environment (not in the centre of a cluster 
or group at 2 = 0). The cosmology employed in the simula- 
tion is Q x = 0.742, On = 0.258, Q h = 0.045 and a 8 = 0.8, 
with comoving box-size = 9 Mpc h _1 and the simulations 
start at z & 499 (p arameters are consistent with WMAP5 
(|Punklev et al.ll2009l )). There are three nested grids located 
in the region where the galaxy halo will form which cor- 
responds to a volume of approximately (2.7 comoving Mpc 
h -1 ) 3 . This has the equivalent resolution of 1024 3 dark mat- 
ter particles which each have a mass of M DM 5.4x10 4 Mq. 
Further levels of refinement, up to a maximum of 15, are ap- 
plied to the grid in this region in order to keep the spatial 
resolution, in the densest region, below ~ 1 pc physical at 
all times. We use a quasi- lagrangian refinement strategy (i.e. 
the gas mass of cells stays roughly constant) in which the 

1 Nut is the Egyptian goddess of the night sky 



next level of AMR is triggered when the baryonic mass in a 
cell reaches 8 x m sp h (where m sp h = 9.4 x 1O 3 M ), or the 
number of dark matter particles in the cell reaches 8. Fig.[T] 
illustrates the grid structure (top row) and corresponding 
gas density with halo and subhaloes indicated with white 
circles (bottom row), showing the set-up of the nested grids 
and the locations where additional levels of AMR are trig- 
gered (see caption for details). Additional technical details, 
including descriptions of the star formation and SN feedback 
prescriptions, are given in the Appendix. 

Fig.[2]shows temperature maps of a region of size 12r v ir 
(top row) and density maps of a region of size 0.2r v ir (bot- 
tom row) centred on the resimulated galaxy at z = 9 for 
the cooling (left column) and feedback (right column) runs. 
Since cooling is present in both simulations, a thin, rapidly 
rotating disc forms, which is continuously fed by cold gas 
that streams along the filaments at ^Mach 5. The disc 
is gravitationally unstable and occasionally fragments into 
star-forming clumps of several parsecs in radius (e.g. bot- 
tom right panel of Fig. [2]) which our resolution allows us 
to resolve along with the filaments and the scale-height of 
the thin disc (this is demonstrated in Powell et al, in prep). 
In our simulation with SN feedback a far-reaching galactic 
wind develops at very high-redshift. Star formation in the 
galaxy, its satellites and haloes embedded in the filaments 
results in many SNe explosions, which are individually re- 
solved. The numerous bubbles overlap, creating an extended 
cavity through which the hot gas escapes the galactic poten- 
tial in the form of a wind. By z = 9, when the dark matter 
halo mass has reached ^ 5 x 10 9 M©, the wind extends to 
~ 6r v ir i-e. ~ 30kpc, filling the virial sphere with hot, diffuse 
gas. 



3 IDENTIFYING THE WIND AND 
FILAMENTS 

In order to quantify the relationship between filamentary ac- 
cretion and galactic winds, we define 5 different gas phases, 
using temperature and density cuts and measure mass inflow 
and outflow rates in each of these phases. 

Fig. [3] shows density-temperature histograms for gas in- 
side the virial radius in the cooling (top) and feedback (bot- 
tom) runs. The region of cold (T < 10 4 K), dense (p>10 
atoms/cc) gas in both phase diagrams can be attributed to 
gas that has condensed into galaxies, either the main galaxy 
or its subhaloes. The cooling is more efficient in the feedback 
run due to the increased metallicity following enrichment 
from SNe. This accounts for the greater mass fraction of gas 
in this region at the lowest temperatures and highest den- 
sities. The edge apparent here is due to the poly trope used 
to prevent artificial fragmentation below the Jean's scale 
(more details in the Appendix). A striking feature, common 
to both the cooling and feedback runs, is the prevalence of 
gas with T ~ 10 4 K. We identify this with filaments (see 
Fig. [4]) defined as containing all gas with T < 2 x 10 4 K, 
0.1 < p < 10 atoms cm -3 . 

Since we are analysing high redshift (z ^ 9) outputs, 
the virial temperature is always below T — 2 x 10 5 K. While 
this temperature threshold has been employed in other stud- 
ies t o successfully sepa r ate shocked-heated gas from cold gas 
(e.g. iKeres et al.ll2005l : lOcvirk et alJl2008h . it does not iso- 
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Figure 2. Top row: Density weighted temperature for a cube of 12r v i r ~ 60 kpc on a side centred on the main halo for the cooling 
(left) and feedback (right) runs at z = 9. Black circles in the top images indicate the position and virial radius of the main halo (the 
largest circle) and its subhaloes. The colour scale has been fixed to 3.0 ^ logT ^5.3; logT = 5.3 is the threshold for gas to be considered 
'hot' (see Table [TJ and the lower limit of logT = 3 has been chosen such that the filaments can be seen clearly since a significant mass 
fraction of the gas is in the range 3.3 ^ logT ^ 4.3. Bottom row: Projected density for a cube of 0.2r v i r on a side centred on the 
main halo, showing the disc approximately face-on in the cooling (left) and feedback (right) runs at z = 9. Note the clumps visible in 
the feedback run (right). 



late all the gas heated to temperatures of the order of T v i r 
and higher in our simulations (as is clear from Fig. [3]). To 
account for the gas which has experienced heating (i.e. has 
T > 2 x 10 4 K and so is not 'cold'), but has not reached 
the 'standard' definition of hot used in the aforementioned 
literature (T > 2 x 10 5 K), we define an additional cate- 
gory of warm gas which has 2 x 10 4 K < T < 2 x 10 5 K. In 
the bottom panel of Fig. [3] (feedback run) we can clearly 
see a cloud of hot (T > T vir « 1.5 x 10 5 K), diffuse (p<10 
atoms/cc) gas which contains a significant fraction of the to- 
tal gas mass. Perhaps surprisingly, we see a similar feature 
in the top panel of Fig. [3] for the cooling run, albeit one that 
makes a less significant contribution to the total gas mass. It 
is predicted that infalling gas will shock on impact with the 
dense ISM if it has not passed through a shock at the virial 
radius (and so was accreted cold) (e.g. Birnboim & Dekel 



2003). The hot phase in the cooling run must be the result 
of shock-heating of this type, since there are no SNe to heat 
the gas in that simulation. However, since the mass in the 
hot phase in the cooling run is tiny compared to that in the 
filaments (see e.g. colours for the hot diffuse gas component 
at z=9 in Fig. [3] (top panel)), and since the cooling time of 
metal poor material at 10 6 K with a density comparable to 
that of the cold filaments (~ 1 at/cc) is around ~ 1 Myr, i.e. 
much larger than the time resolution of the simulation on the 
level of refinement reached at these densities (~ 0.1 Myr), 
we conclude it is primarily the cold diffuse gas that is being 
shock heated. Indeed, if the filamentary material were being 
continuously shocked one would expect to see a prominent 
hot phase with densities equaling or exceeding that of the 
filaments due to compression of the gas as it passes through 
the shock. This scenario is clearly excluded by Fig. [3l which 
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Figure 3. Temperature-density histograms for all gas within the 
virial radius of the main progenitor in the cooling (top) and feed- 
back (bottom) run at z = 9. The colours show the fraction of the 
total gas mass within each bin. Solid horizontal lines show the 
temperature cuts at T = 2 x 10 4 K and T = 2 x 10 5 K and verti- 
cal lines show the density cuts at 0.1 atoms/cc and 10 atoms/cc 
which we have chosen to distinguish the different accretion modes, 
labelled: clumpy ('CL'), filaments ('F'), hot diffuse ('HD'), warm 
diffuse ('WD') and cold diffuse ('CD'). The virial temperature for 
this mass halo (M - 5 x 1O 9 M ) is 1.3 x 10 5 K and 1.6 x 10 5 K 
for the cooling and feedback runs, respectively. 
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Table 1. Summary of the temperature and density thresholds 
used to distinguish different components of the gas inside the 
virial radius of the main progenitor. 



POR (www.vapor.ucar.edu). The results of a p > 10 atoms 
cm -3 cut on the density appears in the top row; the cen- 
tral galaxy and several satellites are clearly denned by the 
clumpy category in both the cooling and feedback runs. The 
gas that meets our definition of the filaments (T < 2 x 10 4 K, 
0.1 ^ p ^ 10 atoms cm -3 ) is illustrated with a render- 
ing of the density in Fig. [4] (middle row). The filamentary 
category is successful in isolating the filaments in both the 
cooling (left) and feedback (right) runs. It appears that the 
filaments are thinner and less well-defined in the feedback 
run, so one already wonders whether this impacts the ac- 
cretion rate. This image also shows an interesting feature 
of the structure of the filaments; they become twisted and 
loop around the central galaxy. This is particularly evident 
in the cooling run (middle row, left) and is linked to the way 
these filaments dissipate their energy as they find their way 
to the centre of the halo. 

Fig. [4] bottom row shows 3D renderings of the tem- 
perature of the diffuse gas (p<10 atoms/cc). It is immedi- 
ately apparent that the virial sphere contains significantly 
more hot gas in the feedback run (right) than the cool- 
ing run (left). Approximately 50 per cent of the virial 
sphere is filled with hot gas in the feedback run, yet only 
5 per cent in the cooling run. This large covering factor 
of hot gas even in haloes below the shock mass threshold 
could explain recent observations of high-redshift galaxies 
in which outflows, rather than cold inflows (as expected for 
cold-accretion dominated haloes) are o bserved (|Steidel et al.l 
bOld iFaucher-Giguere fc Keresl l2010h . The hot gas in the 
cooling run (bottom row, left) surrounds the wound up fila- 
ments and originates from the filaments dissipating their ki- 
netic energy as they stream through the cold, diffuse medium 
on their approach to the centre. We believe this loss of en- 
ergy in a shock in combination with gravitational focusing 
causes the filaments to wrap around the galactic disk. We 
explore the details of this in Powell et al, in prep. 



shows that only a tiny mass fraction of such high density hot 
material exists. We undertake a more detailed study of how 
and where shock-heating occurs in Powell et al, (in prep). 

The categories we have created are summarised in Ta- 
ble [T] and in Fig. [3] we show how the thresholds divide up 
the p — T phase diagrams at z — 9 for the cooling (top) 
and feedback (bottom) runs. Fig. [4] shows visualisations of 
these categories created using 3D visualisation software VA- 



4 THE GALACTIC WIND 

We identify the galactic wind illustrated in Fig. [2] (top, 
right) with both the warm and hot diffuse phases denned 
in the previous section. It is clear from the phase diagrams 
in Fig. [3] that these phases contain significantly more mass 
in the feedback run than in the cooling run and the only 
possible culprit is SNe feedback. We quantify the extent of 
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Figure 4. 3D visualisations (made with vapor) at z = 9 to illustrate the different gas phases in which we measure the different modes 
of accretion and the outflows in the cooling (left) and feedback (right) runs. The cube is centred on the main progenitor halo, extends 
out to r v j r and was extracted on level 15 (equivalent resolution = 40pc). Top row: Rendering of the density for p > 10 atoms/cc showing 
the main galaxy and some satellites. Middle row: Rendering of the density in the range 0.1 < p < 10 atoms/cc for which we define our 
filamentary accretion. Bottom row: Rendering of the temperature for gas with p < 0.1 atoms/cc and T < 2 x 10 4 K representing the 
cold diffuse phase (light-dark blue), and gas with p < 10 atoms/cc and T > 2 x 10 4 K representing the warm and hot diffuse phases 
(yellow-red). Note that gas with T > 2 x 10 5 K is shown in red for the cooling run (left) and orange in the feedback run (right). Colour 
bars show density scale in log atoms/cc (top and middle rows) and the temperature scale in log K (bottom row). The x, y and z unit 
vectors are shown with red, green and blue arrows, respectively. 
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Figure 5. Volume-weighted temperature as a fraction of T v j r in 10 pc radial bins in a sphere of radius 31.8kpc physical (6r v i r at z = 9) 
centred on the halo of the main progenitor, versus redshift for the feedback run. The solid line shows r v j r of the main halo at each 
redshift. The volume- weighted temperature is above T v j r out to ~ 5 — 6r v j r at z < 14 indicating that the wind has escaped the potential 
well of the main halo. 



the galactic wind by examining the time evolution of the 
volume-weighted temperature (represented by the colour- 
coding) as a function of radius in Fig. [5] The temperature 
is calculated in 10 pc physical radial bins within a sphere of 
radius 31.8 kpc physical (6r v ir at z = 9) centred on the main 
progenitor halo at each redshift. By z — 9 the hot gas in the 
feedback run already extends out to ~ 5 — 6r v ir- We note 
that such a wind 'bubble' extent is in f air agreement with 
the es timate derived semi-analytically bv lFurlanetto &; Loebl 
(2003) (their Fig 1, top right panel) using a spherical thin- 
shell approximation. This is rather surprising in light of all 
the simplifying assumptions that go into their modelling. 
For instance, they assume that star formation and feedback 
only occur in galaxies when they cross the cooling threshold 
and are captured by the main progenitor halo. In contrast 
star formation in our simulated galaxies has started earlier 
in most cases and winds have had time to develop prior to 
mergers. The similar sizes of the wind in semi-analytic cal- 
culations arise from the fact that Furlanetto and Loeb's key 
(incorrect) assumption that all the gas has collapsed in a 
dense medium before triggering a starburst and that there- 
fore the wind's spherical thin shell expands into an empty 
halo, actually has consequences very similar to those in the 
scenario which unfolds in our feedback simulation. In the 
latter, most of the halo is also devoid of dense material, but 
this is because this gas is preferentially channelled along 
narrow filaments which continuously fuel star formation, re- 
gardless of whether galaxy mergers are taking place or not. 
Winds from our simulated galaxies expand perpendicularly 
to this filamentary network into the low pressure intrahalo 



medium and depending on their sizes, galactic wind bubbles 
can merge before or after their host haloes merge. 



4.1 Development of the wind 

In order to demonstrate how this far-reaching wind develops 
and to corroborate the picture sketched out in the previous 
section, we examine how the gas inside rw is distributed 
both in temperature and radius in Fig. [6] for the cooling 
(top row) and feedback (bottom row) runs. For each red- 
shift output, we measure the temperature distribution in 
spherical shells of radial width 10 pc physical centred on the 
main progenitor. The colour scale shows the mass fraction 
within each temperature-radius bin and the virial tempera- 
ture at each time is indicated with a horizontal line. As we 
have already seen in Fig. O most of the gas mass in both 
the cooling and feedback runs is concentrated at T ~ 10 4 K 
at all radii. Gas at this temperature spans a large range of 
densities (approximately 0.01 — 100 atoms/cc), however the 
main mass contribution is from the filaments (see Section [3] 
for more details). The most striking difference between the 
runs is the prevalence of hot gas (T > T v i r ) in the feed- 
back run (bottom row) at all radii, in stark contrast with 
the cooling run (top row); this gas constitutes the galactic 
wind. 

There are many 'spikes' in temperature in both the cool- 
ing and feedback runs, which correspond to the gas in sub- 
haloes or other haloes beyond the virial radius of the main 
halo. We demonstrate this explicitly for subhaloes within 
r v ir of the main halo, by overplotting squares at T = 10 4 K, 
to indicate the position of the centres of the subhaloes with 
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Figure 6. Time evolution of the temperature-radius histograms for all gas out to 6r v i r in the cooling (top row) and feedback (bottom 
row) runs. The colours show the logarithmic fraction of the total gas mass within each bin. The position of the main halo (by definition 
at radius= 0) and its subhaloes are overplotted with light green squares at T = 10 4 K, with bars to indicate 0.1r v i r for each. A solid 
black horizontal line marks the virial temperature, specified in the top right corner of each plot. 




Figure 7. Time evolution of the temperature-radius histograms for all gas out to 6r v i r in the feedback run. The colours show the average 
metallicity in each bin. The position of the main halo (by definition at radius= 0) and its subhaloes are overplotted with light green 
squares at T = 10 4 K, with bars to indicate 0.1r v i r for each. 



respect to the centre of the main halo (which is at r = 0, by 
design). The horizontal error bars indicate 0.1r v ir for each. 
The low temperature end of the spike indicates where gas 
has cooled and condensed within the subhaloes. 

The high temperature ends of these spikes are of partic- 
ular interest in the context of the galactic wind. In Fig. [71 we 
show the T — r diagram from Fig. [6] but now the colour scale 
indicates the mass-weighted metallicity. The high mass- 
weighted metallicity we measure in the high-temperature 
spikes corresponding to (sub) haloes shows that this hot gas 
results from SNe exploding in the galaxies hosted by these 
haloes. By examining the time sequence of plots for the feed- 
back run (Fig. [6] bottom row and Fig. [7]), we can see that 



these plumes of hot, metal-enriched gas (seen particularly 
clearly at z = 12.3, second panel of bottom row of Fig. [6] 
and Fig. [7} gradually build up into a continuous distribu- 
tion of hot gas at all radii, forming the galactic wind. 

In Fig. [8] we show one cell thick (~ 0.05r v ir) slices 
through the gas density in a region of size 12r v ir x 12r v ir, 
centred on the main halo. Here we can clearly see the low 
density cavities left behind as the blastwaves from the SNe 
propagate outwards. At higher redshifts (towards the left) 
these 'bubbles' originate only from the main galaxy, yet as 
we move to lower redshifts (towards the right), cavities can 
also be seen forming further out in the filaments as galax- 
ies in haloes embedded therein also host SNe explosions. By 
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z — 9 (the last panel) we can see the bubbles are overlapping 
merging into one outflow with its epicentre located at the 
main galaxy. 

In previous studies where galactic winds were in- 
voked, often the wind could not escape the galaxy (e.g. 



Ocvirk et al. 20081). unless i t was put in by han d 



Springel fe Hernauistl 120031 : lOppenheimer fc Pavel 12008 
We propose it is this mechanism of overlapping bubbles 
from multiple (sub) haloes which (i) significantly increases 
the reach of the wind (ii) more efficiently pushes material 
outside the main halo. Indeed, if one considers that the to- 
tal amount of kinetic energy imparted to the gas by all SNe 
(those exploding in the main halo and its subhaloes) is fixed, 
but that one is allowed to distribute it at different points 
along the path of the outflow, rather than release it all in 
the central galaxy, then more energy is available to drive ma- 
terial out of the main halo's potential well. The reasons are 
twofold: (i) the SN driven winds of satellite galaxies have 
to climb out of shallower potential wells (ii) the (smaller) 
amount of gas swept from the central object (SNe exploding 
in subhaloes will not contribute to the mass loading) will see 
its energy replenished on the way outwards by interacting 
with the subhalo winds. Note that as a result of this effect, 
galactic winds might carry more energy and less mass out 
of the halo. 

High temperature spikes coincident with a (sub) halo are 
also occasionally evident in the T — r diagrams for the cool- 
ing run, however these are of a similar colour but much thin- 
ner than in the feedback run and therefore contain a much 
smaller fraction of the gas mass. This hot gas is generated in 
the shocking of the diffuse circum-galactic material of satel- 
lite galaxies (see Powell et al, in prep for a detailed discussion 
of shock- heating in the cooling run). 

Another feature of the T — r diagrams (Figs. [6]and [7} is 
a haze of gas (i.e. gas not associated with one of the spikes) 
with T < 10 4 K which persists at all radii. This is cold, 
metal poor, diffuse gas that fills the space between the fila- 
ments. In the feedback run (bottom row of Fig. [6} this gas 
only becomes significant beyond about 3r v ir (~ 15kpc) at 
z = 9, whereas it is present in significant quantities within 
r v ir in the cooling run (top row). Furthermore if we follow 
the plots for the feedback run (bottom row of fig. [6]) from 
high to low redshift (left to right) we can see this haze of cold 
diffuse gas recedes over time. The hot wind from the SNe is 
sweeping away the cold diffuse gas which surrounds the fila- 
ments, reducing the resistance to future outflows from SNe. 
We substantiate this hypothesis in Section [5] with detailed 
measurements of the outflow. 

4.2 Wind properties 

In Fig. [9] we present the properties of the galactic wind as a 
function of radius (within r v ir of the main halo) at z = 9. We 
show the mass- weighted velocity (top left), the metallicity 
profile (top right), the net mass outflow rate (bottom left) 
and the net mass outflow rate of metals (bottom right). The 
solid red line shows these quantities for the warm and hot 
gas combined, which is then split into warm (orange dashed 
line) and hot (red dashed line) components (see Table [1] for 
definitions of 'hot' and 'warm'). 

The velocities we measure (top left p anel of Fig. Q com- 
pare well to the semi-analytic model of iFurlanetto &; Loebl 



(2003), both for launching (their equation (11) for v w ) and 
average expansion velocities (equation (12) and Fig. 2 top 
right panel for v exp ) of the wind. Agreement on v exp is not 
too surprising because the wind bubble in our simulation is 
able to avoid zones of high pressure (filaments have a neg- 
ligible halo volume filling fraction), and therefore it mostly 
sweeps gas with the average density (and pressure) of the 
IGM at the redshift of interest. This makes Furlanetto and 
Loeb's approximation of a spherically symmetric expansion 
in a homogeneous background of average density a good one, 
even though it is not very consistent with the spherically 
symmetric accretion of dense gas which should naturally oc- 
cur in their model if they did not wait until all this gas 
has piled up in the centre of the halo before triggering their 
wind. 

More remarkable is the agreement with v w , consider- 
ing that the fraction of the total energy of SNe available to 
power the wind in their model (/esc) is calibrated on the ide- 
alized simulations of high redshift starbursts by iMori et"all 
(2002) who quote a value of about 25%. Indeed, using the 
same estimator as Mori and collaborators (see Fig. 12 of 
their paper), i.e. the ratio of kinetic energy flux at the virial 
radius to the total supernova luminosity, we find that / esc is 
around ~ 2.5%. This remains the case even when we vary 
the radius at which we measure the kinetic flux (down to 0.1 
r v ir) or the time delay (10 — 100 Myr) between the measure- 
ment of the kinetic flux at the virial radius and the super- 
nova luminosity to account for the amount of time it takes 
the central wind to reach the growing virial radius in our 
non-idealized feedback simulation. This factor 10 lower for 
the efficiency to power winds in our simulation compared to 
that measured by Mori and collaborators can be attributed 
to (at least) two diff erences: (i) the central density of the gas 
in lMori et all (|2QQ2h is severely underestimated (they use 10 
at / cc which is at the bottom end of our definition of clumpy 
gas and more representative of the density of filaments at 
these redshifts) (ii) the star formation rate (and hence SNe 
explosions) is spread out in time (and space) in our sim- 
ulation rather than instantaneous. These two effects, and 
especially the former, mean that radiative cooling affects 
our supernova remnants to a much greater level, sapping 
their energy much more efficiently. Finally, we note that the 
value of /esc in our simulation does not seem to be a strong 
function of the mass of the halo, as the progenitor of our 
5 x 1O 9 M0 halo at z = 9 reaches a similar mass as that 
of the halo used by Mori and collaborators (2 x 1O 8 M0) at 
z ~ 13 without /esc changing much. Thus, the agreement 
with IFurlanetto &; Loebl (|20Q3h for v w (the wind launching 
velocity) is recovered because they assume that the the ef- 
ficiency of the wind mass loading, / sw , is related to the star 
formation efficiency on galactic scales, /*, by / sw = 2/*. Al- 
though the values for the /* that we measure are close to 
the the value of 10 % that these authors adopt (see section [6] 
for details), the value we measure for / sw , is approximately 
a factor 10 lower than theirs (see section [5] for the exact 
value). Since we adopt an initial mass function (IMF) very 
similar to theirs (Salpeter with similar mass cut-offs), the 
order of magnitude discrepancies in / esc and / sw cancel out 
in their formula (11) for v w and we recover a value close to 
their 300 km/s, as shown in the top left panel of our Fig [9] 

The metallicity profile of the wind out to 6r v ir in Fig. [9] 
(top right panel) shows that the metallicity ranges from 
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Figure 8. The gas density in a z-slice of thickness 0.05r v i r and size (12r v i r ) 2 centred on the main progenitor for the feedback run at 
(from left to right) z = 15.7, 12.3, 10.8 and 9. Arrows show the momentum vectors (pv). The low density cavities are created as the SNe 
blastwaves propagate outwards and eventually overlap producing a far-reaching galactic wind. 



0.5^0 to O.IZq at z = 9. In terms of mass, 75 per cent of 
the gas with T > 2 x 10 4 K within r V i r has Z > O.IZq, with 
quite narrow scatter since only 10 per cent of this gas has 
Z > O.5Z0. Whilst we cannot currently compare measure- 
ments of the metallicity of a high-redshift galaxy directly 
with observations, we can draw parallels with observations 
of the metallicity of winds in local galaxies. Determining the 
wind metallicity using X-ray observations of the hot gas is 
difficult due to va rious uncertainties and degeneracies (see 
IVeilleux et alj|2005l for a dis cussion). Bearing these caveats 
in mind jMartin et al.l (|2002h find that the metallicity of the 
wind in NGC 1569, a local dwarf starburst galaxy, must 
be greater than 0.25Zq and that the wind seems to carry 
most of the metals produced in the current starburst. Our 
measurements of the metallicity of the hot diffuse gas are 
broadly consistent with these observations, suggesting that 
metal enrichment away from the main site of star formation 
is plausible at such a level. 

If we extrapolate our results using Press -Schechter the- 
ory to calculate the average number density of 5 x 10 9 Mq 
haloes at z = 9, we find a value of « 0.15 Mpc -3 comov- 
ing for the cosmological model we use. This translates into a 
mean inter-halo distance of about 190 physical kpc, so that a 
small but non negligible volume fraction (0.4%) of the IGM 
can be filled by metal rich hot gas wind bubbles blown by 
haloes in thi s mass range at those re dshifts. Assuming that 
the result of iFurlanetto &; Loebl (|2003h holds in this regime, 
such that the extent of the wind does not depend sensitively 
on the virial mass of the galaxy host halo at a given redshift, 
then the volume fraction filled by hot gas will increase if we 
start including lower mass haloes. Based on the measure- 
ment of the wind metallicity in our simulation (see Fig. [9] 
top right and bottom right panels), it seems that galactic 
winds may indeed be very efficient at polluting gas in the 
the IGM at very high redshifts up to a level comparable to 
that observed in the Ly-a forest at z ~ 2 (jBouche et al.l 
l2007h . 



summarised in Table [TJ We now measure the rate of mass in- 
flow and outflow in each component in order to quantify how 
they contribute to the overall accretion/outflow of material 
onto/from the central galaxy. In particular, we focus on es- 
tablishing the impact (if any) of the galactic wind (which 
we identify with the warm and hot diffuse phases) on these 
mass fluxes by comparing the feedback and cooling runs. 

We measure the flux in spherical shells centred on the 
centre of the dark matter halo of the main galaxy, which is 
taken to be the position of the densest dark matter particle. 
The shells are of width 100 pc physical and extend out to 
the halo's virial radius. In order to compute gas velocities 
relative to the galaxy, we subtract the peculiar velocity of 
the dark matter halo from all gas velocities. We compute 
the mass flux, rafl ux , for each grid cell in a shell using, 



■ Pcell ^radial G^E 



(1) 



where p C eii, ^radial and dx are the cell's gas density, radial 
gas velocity, and physical size respectively. 

By definition, gas with negative mfl ux (i.e. negative 
^radial) is flowing inwards and gas with positive rafl ux (i.e. 
positive Radial) is flowing outwards. We compute the total 
inflow, F-, and total outflow, F+, separately for each shell 
as follows, 



£ra flux , + 2 



(2) 



where r s heii is the distance between the halo centre and 
the shell midpoint and the +/— subscript on rafl ux denotes 
whether it is positive or negative. This calculation is re- 
peated for each of the gas components listed in Table [1] the 
only difference being that the terms rafl ux ,+ and mflux,- only 
include contributions from cells that meet the temperature 
and density criteria for the particular component. The flux 
we compute is, therefore, an instantaneous flux at the time 
of the output we are analysing. 



5 MEASURING GAS ACCRETION AND 
OUTFLOW 

In Section [3l we defined 5 gaseous components (clumpy, fila- 
ments, cold diffuse, warm diffuse and hot diffuse) which are 



5.1 Inflow 

Fig. [10] shows the gas inflow versus radius (up to the virial 
radius) for the cooling (first column) and feedback (second 
column) runs. The mass flow is measured for each gas phase 
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Figure 9. Mass- weighted radial velocity profile averaged (top left), metallicity profile (top right), net mass outflow rate (bottom left) 
and net mass outflow of metals (bottom right) in 100 pc (physical) spherical shells out to 6r v i r for the feedback run at z = 9. Solid red 
line shows quantities warm + hot gas, which is then split into warm (orange dashed line) and hot (red dashed line). 



in spherical shells of width 100 pc physical, centred on the 
main progenitor halo, out to its virial radius. 

The total mass inflow rate (solid black line) is ap- 
proximately constant with radius, reaching values above 
IMoyr -1 by z < 12, in both runs and is dominated by the 
contribution from the filaments (blue dotted line). This level 
of inflow compares well with the estimate of the baryonic ac- 
cretion rate at the virial radius derived by simply painting 
baryons onto dark matter using extended Press- Schechter 
theory (|Neistein et alJl2006h . For ACDM, this yields: 

(jDekel et al.l l2009h . In this formula, fh is the universal 
baryon fraction which is 0.17 in our simulation, z — 9 and 
Mvir « 5 x 1O 9 M , so we get M 2.8M yr _1 , which is 
within 50 percent of our measured values in the cooling and 
feedback runs (cf. top two rows of Fig. Hip . At z — 15.6, 
the virial mass of our dark matter halo has dropped by a 
factor 60, to M v i r ~ 8 x 1O 7 M . Using the same formula 
we therefore predict an accretion rate of M — O.O73M0yr _1 
at this redshift, which is lower, but of the same order of 
magnitude to that which we measure (M O.SMoyr- 1 ) m 
the simulation. This confirms that our resimulated halo is 
a fairly 'typical' halo and that gas, even in high resolution 
simulations is being accreted proportionally to dark matter, 
at least at the virial radius. 

While the similarity of the gas accretion rate to the pre- 
diction for the dark matter accretion rate at the virial radius 



is not very surprising, the near constancy of the filamentary 
gas accretion rate all the way down to the disc is more re- 
markable. This is caused by the radial velocity of the gas 
in the filaments remaining constant around ~ 50km/s, in- 
dependently of radius. Indeed, whereas the dark matter fila- 
ments are dynamically heated when they penetrate the virial 
radius, the gaseous filaments radiate away the extra grav- 
itational energy (recall their temperature is stable around 
~ 10 4 i^) to continue their journey supersonically (at about 
Mach 5) deep into the centre of the dark matter halo, all the 
way to the central disc. 

On the other hand, the mass inflow rate in the cold dif- 
fuse component (blue dashed lines Fig fTU]) increases with ra- 
dius, until it becomes comparable with that in the filaments 
at the virial radius. This trend with radius is a geometrical 
effect; as the shell surface area decreases with decreasing ra- 
dius, more and more diffuse gas is compressed and joins the 
filamentary phase, and thus the inflow of diffuse cold gas per 
unit time is reduced as the radius shrinks. As the filaments 
have a roughly constant cross-sectional area (as illustrated 
in the middle row of Fig. [4|), no such effect occurs and thus 
they rapidly become the dominant supply of gas in the inner 
regions of the halo. The accretion rate is very similar for the 
feedback and cooling runs and so it seems that the geometry 
of the supply of cold gas to the central object means that 
the SNe feedback is unable to significantly diminish it. The 
main source of such gas is the dense, collimated filaments, 
which present more pressure resistance to the galactic wind 
than the more symmetrically distributed diffuse material. 

The impact of SNe feedback on accretion is also con- 
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Figure 10. Gas inflow rates averaged in 100 pc (physical) spherical shells out to r v - 1Y for the cooling (first column) and feedback (second 
column) runs. Black solid lines show the mass inflow rates of all gas, which is then split into dense (blue dot-dash line), filamentary (blue 
dotted line), cold diffuse (blue dashed line), hot diffuse (red dashed line) and warm diffuse (orange dashed line) categories as summarized 
in Table 1. The solid red line shows the warm diffuse and hot diffuse components combined. 
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Figure 11. Gas outflow rates averaged in 100 pc (physical) spherical shells out to r v - ir for the cooling (first column) and feedback (second 
column) runs. Black solid lines show the mass outflow rates of all gas, which is then split into dense (blue dot-dash line), filamentary (blue 
dotted line), cold diffuse (blue dashed line), hot diffuse (red dashed line) and warm diffuse (orange dashed line) categories as summarized 
in Table 1. The solid red line shows the warm diffuse and hot diffuse components combined. 
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Figure 12. Left: Star formation rate averaged over 10 Myr intervals versus redshift for the feedback run (dashed red line) and cooling 
run (solid blue line). This is calculated using the stars within the virial radius of the main progenitor at each redshift and indicates the 
in-situ star formation rate for the main progenitor itself. Right: Star formation of main progenitor only (averaged over 10 Myr intervals) 
for the feedback run divided by that for the cooling run, as a function of redshift. The horizontal line indicates when the SFR in the 
feedback run equals that in the cooling run. 



sidered in Ivan de Voort etHl (|2010h . in which the authors 
examine accretion onto a large sample of galaxies and their 
haloes in a suite of cosmological simulations with different 
physics. We note that their recipe for SNe feedback involves 
adding in a galactic wind 'by hand' (i.e. by applying a ve- 
locity kick to particles near an exploding star particle) and 
so their technical approach is very different to ours. A di- 
rect comparison to our results is challenging, not least be- 
cause we are considering a very different halo mass and 
redshift range (~ 5 x 1O 9 M0 halo at z > 9, compared to 
haloes with M ^ 1O 1O M0 at z < 5.5). In the discussion 
that follows, we will compare to their results for haloes with 
M vir <2 x lO n M0 at z — 2, as this is the closest match to 
the halo in our study. Ivan de Voort etaD (|2010h find that 
including SNe feedback reduces the total smooth accretion 
rate by a factor of « 2.5 for haloes in this range and by a 
factor of ~ 10 for the central galaxies hosted by these haloes. 
This suggests that accretion from the IGM has been consid- 
erably impeded by SNe feedback; this is in stark contrast to 
our results. 

A clue to the possible origin of this difference can be 
found in their Fig. 7 (last row) where the lowest mass halo 
illustrated (which has 10 11-5 M©) is embedded within a fil- 
ament. This means there is no lower density material for a 
SNe-driven wind to preferentially flow into. This is very dif- 
ferent to our lower mass, but higher redshift halo, which is 
fed by 3 filaments that are thin compared to its virial radius. 
It is interesting to note that for their 1O 12 M0 halo, which 
is fed by 3 filaments (see their Fig. 7, middle row) the total 
accretion rate onto the halo is not altered by SNe feedback. 
We would therefore argue that the geometry of the accretion 
must be considered, not just the halo mass (although the 
two are interlinked) in order to accurately predict /explain 
the effects of galactic winds on gas accretion. 

The spikes in the inflow in the clumpy component (blue 
dot-dash lines in Fig |10[) correspond to satellite galaxies, or 
other over-dense clumps of gas. Note that there is always 



a spike in the clumpy inflow rate at the centre (at r = 0); 
this is due to the presence of a rotating gas disc, which also 
gives rise to similar central spikes in the corresponding out- 
flow plots (Fig. 1111 Indeed, since gas discs are generally not 
axisymmetric, as different rotating elements of gas of the 
same disc pass through a spherical shell they can contribute 
either negatively or positively to the radial velocity compo- 
nent of the shell, therefore leading the gas in the disc to be 
detected both as inflowing and outflowing. 

There is both warm inflow (orange dashed line) and to 
a lesser extent, hot (red dashed line) inflow, in the cooling 
run (Fig. \10\ first column). While this is around 3 orders 
of magnitude less than the contribution from the filaments 
it indicates the possibility for gas to be heated in low mass 
haloes even without feedback. Fig. U revealed that the origin 
of this heated material is a shock front, where cold diffuse 
material is heated as it impacts the denser filamentary ma- 
terial surrounding the disc (this, and the properties of gas 
inflow in general, is discussed in more detail in Powell et al, 
in prep). There is a similar rate of mass inflow in the hot 
and warm components in the feedback run, but it is dwarfed 
by the rate of outflow. 



5.2 Outflow 

Fig. II II shows the gas outflow versus radius (up to the virial 
radius) for the cooling (first column) and feedback (second 
column) runs, respectively. The rate of mass outflow in the 
feedback run is dominated by the contribution from the hot 
component. Despite velocities reaching several 100 km/s (see 
Fig. [9] left) for this component, the mass outflow rate is only 
around 10 — 30 percent of the total inflow. This is simply 
because the hot gas has very low densities; Fig. [3] (bottom 
panel) shows that the greatest mass fraction of hot gas oc- 
curs at densities of ~ 10 ~ 3 atoms/cc. 

The cold diffuse outflow in the feedback run is a factor 
of 10 higher than that in the cooling run (see particularly 
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z = 9) , reinforcing our interpretation of Figs [6] and [71 that 
hot outflows from SNe sweep away some of the cold diffuse 
gas falling onto the central galaxy. This fraction of swept- 
up material is about 10 percent of the amount of cold dif- 
fuse gas which flows in, suggesting that the material being 
swept out is the low density (<~ 10 ~ 3 atoms/cc) cold (10 3 K 
<T < 10 4 K) tail of the distribution seen in the p — T phase 
diagram for the cooling run (Fig. [3] top panel). Indeed this 
gas is clearly missing in the feedback run (compare the sec- 
tions labelled 'CD' in Fig. [3) for the cooling (top) and feed- 
back (bottom) runs). As with the filamentary component, 
the impact of feedback on the cold diffuse accretion rate is 
surprisingly small even though a high velocity galactic wind 
has been successfully launched. 

The large spikes in the mass outflow rate of the clumpy 
component (blue dot dash line) in both runs are due to the 
presence of rotating gas discs in satellite galaxies. The out- 
flow seen in the filamentary component (blue dotted line) is 
a result of the 'looping' of the filaments (particularly pro- 
nounced in the central regions of the halo in the cooling run) 
and, to a lesser extent, material associated with the outer 
region of satellites which is not distinguished from the fila- 
ment in which the satellite is embedded due to the smooth 
density gradient. 



6 CORRELATING ACCRETION, OUTFLOWS 
AND STAR FORMATION 

One of the main reasons why the cold-mode of accretion has 
received so much attention, is because of its implications for 
the star formation histories of galaxies. In the traditional 
picture, star formation is delayed as the gas is first shock 
heated to the virial temperature of the dark matter halo 
and requires time to cool and condense isotropically onto 
the central galaxy. If, on the other hand, gas comes in cold 
and fast along dense filaments, star formation can commence 
earlier, be much more rapid and localized in space. In order 
to quench such starburst episodes, the supply of fuel for star 
formation has to be shut off by e.g. SNe feedback. However, 
since galactic winds are expected to take the path of lowest 
pressure to propagate in the halo, the filaments could be 
resistant to destruction and the starburst prolonged. There- 
fore, in this section we investigate how star formation rates 
correlate with net mass inflow rates, comparing the cooling 
and feedback runs to quantify the impact of SNe. 

To examine the star formation rate (SFR) of the main 
progenitor in isolation we analyse the stars within the virial 
radius at each available output and compute the SFR in a 
small time interval only (10 Myrs), thereby excluding stars 
formed outside the main progenitor and accreted during 
mergers (left panel of Fig P"2]) . The right panel of Fig. [12] 
shows the ratio of these in-situ star formation histories in 
the feedback and cooling runs. Whilst the SFR in the cool- 
ing run is a factor 2 larger than that in the feedback run for 
z > 12, the situation is reversed at lower redshifts (z < 10) 
where the SFR of the feedback run exceeds that in the cool- 
ing run by up to a factor 10. This is intriguing as it suggests 
that SNe are a source of positive feedback, i.e. they enhance 
star formation rather than inhibit it. Although in order to 
confirm the origin of this positive feedback effect it would 
be necessary to perform another nut simulation with su- 
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Figure 13. Top: The Kennicutt-Schmidt star formation effi- 
ciency versus redshift for the cooling (blue solid line) and feed- 
back (red dashed line) runs. These quantities are computed for 
the main progenitor, within a sphere of radius 0.1r v i r - Bottom: 
Ratio of net clumpy inflow to net filamentary inflow versus red- 
shift for the cooling (blue solid line) and feedback (red dashed 
line) runs. The value of the inflow rate at each z is obtained by 
finding the net inflow per radial bin and averaging over the radial 
bins within r v i r (full data with radial dependence is presented in 
Figs. fTOl and fTTT). 



pernova feedback turned on but metal enrichment switched 
off, we attribute most of it to the extra cooling provided by 
metal lines in the feedback run rather than, for example, 
compression of the ISM by blast waves. Indeed it is obvious 
from Fig [3] that such an extra channel for cooling enables 
a larger mass fraction of cold gas to reach the star forming 
density threshold in the feedback simulation (bottom panel) . 

To assess how these differences affect the global effi- 
ciency of star formation, we estimate the 'effective' star for- 
mation efficiency, e s f, in both runs by writing a Kennicutt- 
Schmidt law for the entire central galaxy which gives, 

£sf = SFR^- (4) 
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where m gas is the mass of gas measured inside a sphere cen- 
tred on the galaxy, with radius 0.1r v i r and we take td yn ~ £ff 
as the dynamical time. Here ts is the usual free-fall time, 
i.e. the time for the mass enclosed in the sphere to collapse 
under its own gravity i.e., 

The dynamical time varies by a factor of around 2, over the 
redshift range examined (z « 9 — 15), from & 5 Myr at high 
redshift to 10 Myr at low z. 

Fig. [13] (top) shows the efficiency, e s f, as defined in 
equation [4] for the cooling (solid line) and feedback (dashed 
line) runs. Note that this effective (global) star forma- 
tion efficiency is, on average, an order of magnitude higher 
than both the typical star formation efficiency of 0.01 — 
0.02 determined from observations of the local Universe 
(|Krumholz fc TaiJl2007h and the local efficiency parameter 
in RAMSES, which we set to 0.01 in both cooling and feedback 
simulations (for details of the star formation implementa- 
tion see the Appendix). Such low efficiencies of this order of 
magnitude are only measured in the cooling run for z ~ 9. 

This effect can be understood as follows. Equation 2] 
represents an averaging of the star formation properties of 
the whole region we are considering, whereas the efficiency 
parameter in the simulation can only control the star forma- 
tion rate in individual cells which have exceeded the required 
density threshold. If star formation occurs in a clumpy 
medium where many regions exceed the density threshold 
for star formation, the overall efficiency of the object will 
be larger than if the medium is more homogeneous. Higher 
densities in a clumpy medium also give rise to shorter star 
formation time-scales on a cell by cell basis. It seems appeal- 
ing to conclude that, at very high-redshift, star formation on 
galactic scales is more efficient, not because it is more effi- 
cient on local cloud core scales, but because a larger fraction 
of the ISM gas mass is locked in dense clumps (clumps are 
evident by eye in the gas disc in the feedback run: see Fig. [2] 
bottom right), than at low (z ~ 0) redshift. 

This view receives some suppo rt from observations of 
clum py galaxies at z > 1 (e.g. lEhnegreen &; Elmegreenl 
2005). These clump masses are estimated to lie around 
M rsj 1O 8 M0, which are consistent with the fin dings of high 
resolu tion hydrodynamical simu lations. Indeed I Agertz et al.l 
( 2009|) and lCeverino et all (|2010h report massive clumps of 
order 1O 8 M at z = 2 — 3 The clumps we see in our sim- 
ulated disc at z = 9 are less massive (~ 1O 5_7 M ), but 
we believe they are the high redshift counterparts of the 
observed (and simulated) clumpy galaxies. We defer a de- 
tailed study of these stellar clumps to a companion paper 
(Slyz et al, in prep.), but we note that their typical masses 
(Mdump ro 1O 6 M ) are consistent with the Jeans mass val- 
ues one expects from the gravitational fragmentation of the 
gaseous disc. 

In order to understand how the cold gas components 
fuel star formation and how they yield the high global star 
formation efficiencies measured, we evaluate how the SFR 
correlates with the net mass flow averaged over all radii. 
In Fig. [14] (left panel), we examine the ratio of the SFR to 
the cold gas inflow rate versus redshift for the cooling (blue 
solid line) and feedback (red dashed line) runs. In the inter- 



val 9 ^ z ^ 10 the cold gas supply far exceeds the SFR in 
the cooling run. In the feedback run during the same time 
period, the dashed line oscillates around 1, showing that a 
amount of gas comparable to the amount of incoming ma- 
terial is converted into stars. However, the high efficiency of 
star formation that we measure on the scale of the galaxy 
(Fig I13j) seems uncorrelated with the epoch of high inflow 
rates, whereas the global star formation rate is, as shown 
by the lines on the right panel of Fig. [14] Once again, we 
interpret this as evidence that the global efficiency is more 
sensitive to the local microphysics than the global proper- 
ties of mass accretion. Indeed, the epochs of high global 
efficiency correlate quite well with the times at which the 
radially averaged clumpy component of cold gas accretion is 
high (see Fig. [13] bottom), i.e. mergers/strong interactions 
with neighbouring galaxies take place. These are thought to 
trigger violent instabilities which will compress the ISM (or 
increase turbulence and therefore increase fragmentation) 
and allow a larger fraction of it to reach the threshold for 
star formation, resulting in higher global star formation ef- 
ficiencies and a merger induced starburst. Note that when 
the ratio of clumpy to filamentary inflow reaches its peak, 
we also see a corresponding peak in the SFR (left panel, 
Fig. [T2|) . Compared to these violent events the impact of 
smooth filamentary accretion on the global star formation 
efficiency is a second order effect, even though it maintains it 
at a higher level (between 5 — 10%) than the « 1 % expected 
from our choice of efficiency parameter on small scales, which 
should remain unchanged by diffuse accretion alone. 

Finally, Fig. [14] (right panel) shows that the hot outflow 
rate is on average one order of magnitude below the SFR. 
This is contra ry to observations of l ower redshift galaxies in 
iMartinl {l999), and lHeckman et al.1 {2000) who estimate the 
mass-loss rate to be comparable to the SFR. This means 
that the hot gas outflow rate is also about a factor of 10 
lower than the rate of mass entrained by SNe ejecta in our 
modelling (see the Appendix for details of SNe feedback), 
suggesting that a significant proportion of this gas mass does 
reform stars on very short timescales. Whether this is caused 
by a lack of implementation of physical processes known to 
occur in the interstellar medium (e.g. radiative transfer, stel- 
lar winds) or this is an intrinsic problem of (high-redshift) 
SNe driven winds remains to be established. 

In light of this study, we conclude that SN explosions 
alone are very inefficient at driving massive galactic winds 
from typical L* galaxies at high redshift. We stress that 
this finding is resolution-dependent; at lower resolutions 
the Sedov- Taylor solution necessarily represents multiple SN 
remnants and a more massive wind may be driven. Since 
such a solution has only been proven to accurately describe 
an individual SN remnant, we are confident that the low 
mass wind we produce is a robust result. However, there are 
additional processes that could affect the mass-loading of 
the wind which are not included in our simulations, such as 
stellar winds, ionizing photon heating, cosmic rays and mag- 
netic fields. Furthermore, as in all cosmological simulations, 
we use a (simple) subgrid model for star formation and have 
assumed an IMF, both calibrated on local Universe observa- 
tions. As the rate and energy of SNe are dependent on these 
assumptions so, in turn, should be the mass outflow rate in 
the wind. 
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Figure 14. Left: Star formation of main progenitor only (averaged over lOMyr) divided by net cold gas mass inflow rate versus redshift 
for the cooling (blue solid line) and feedback (red dashed line) runs. The horizontal line indicates when the cold gas inflow rate equals the 
SFR. Right: Net mass outflow rate of warm + hot gas divided by the net mass inflow rate of cold gas (solid line) or SFR (dashed line) 
for the feedback run. The value of the inflow or outflow rate at each z is obtained by finding the net flow per radial bin and averaging 
over the radial bins within r v j r (full data with radial dependence is presented in Figs. [T0l and lll[) . 



7 CONCLUSIONS 

We have undertaken a detailed analysis of a suite of 
subparsec resolution simulations (the nut simulations) 
with the aim of understanding the mechanism via which a 
far-reaching galactic wind can arise in a protogalaxy at high 
redshift (z ~ 9) and how this impacts the surrounding IGM 
and the galaxy's evolution. In particular we examine the 
relationship between cold-mode, filamentary accretion and 
star formation and explore the impact that the introduction 
of SNe feedback and the resulting wind has on these 
processes. 

Our main findings are: 

(i) A far-reaching SNe driven wind (illustrated in Fig. [21 
top right) does develop; this is particularly significant since 
we do not directly model the wind in our feedback scheme, 
but instead we resolve the Sedov blastwaves of individual 
SNe and the wind arises naturally. We find that SNe explo- 
sions occur not only in the main progenitor, but also in sub- 
haloes and neighbouring dark matter haloes (demonstrated 
in the temperature-radius histograms, Fig. [6}. These cre- 
ate overlapping bubbles which are critical for creating the 
far-reaching wind and for facilitating the escape of the gas 
from the halo potential. Since subhaloes are key, we specu- 
late that similar winds can also be produced by higher mass 
haloes at lower redshifts, and we will certainly investigate 
this issue in the very near future. The wind is enriched with 
metals to the order of 0.1 Zq (see the metallicity profile in 
FiglU top right). While one must be careful when compar- 
ing high-redshift simulations results to local observations, 
we conclude that it seems feasible for such a wind to enrich 
the IGM to the level required by observations at z « 2 . On 
the other hand, we find the ejection of mass by the wind 
is very inefficient: the outflow rate in the hot phase is only 
10 — 30 per cent of the total mass inflow rate. We also find 



that the mass outflow rate is about an order of magnitude 
lowe r than the SFR, in disagreement wi th local observations 
fe.g |Martinlll999l ; iHeckman et al.ll200d ), perhaps indicating 
an important difference between low and high redshift SNe 
driven winds or a missing physical process in the simulation 
which drastically alters the mass loading factor of the wind 
(stellar winds, radiative transfer). We plan to address this 
issue in a forthcoming paper. 

(ii) The total gas accretion rate at the virial radius of the 
main progenitor (see Fig. 1111 top two rows) ranges between 
1 — lOM yr _1 depending on redshift, and is in fair agree- 
ment with analytical predictions based on extended Press- 
Schechter theory, indicating that the gas accretes like dark 
matter at the virial radius (note that the accretion is cold). 
However, in contrast to the accretion of dark matter, this 
rate persists for the gas filaments right down to the central 
object and thus these filaments are the dominant supply of 
cold gas to the central disc (compared to clumpy or diffuse, 
spherical cold accretion). While the filaments appear more 
perturbed in the inner halo in our 3D visualisations of the 
feedback run than in those of the run without SNe feedback 
(Fig. 2] second row), we measure similar mass accretion rates 
in both runs, suggesting that the majority of the gas survives 
its journey to the central object intact and gets replenished 
by cold gas coming in from larger distances. So, despite the 
development of a galactic wind, the cold, filamentary accre- 
tion is not significantly altered. This suggests that, at least 
at z ~ 9, SNe driven winds cannot reduce the amount of 
cold, accretion onto dwarf-mass protogalaxies and therefore 
cannot significantly reduce star formation at this epoch, as 
is commonly assumed. 

(iii) At the lower end of the redshift range studied (z < 
12), the SFR in the feedback run is more often than not 
greater than that in the cooling run (see Fig. 1141 left for the 
ratio of SFRs versus redshift). An extreme example is at 
z = 9 when the SFR in the feedback run is w 3M yr _1 , yet 
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only O.3M yr _1 (the minimum SFR over the redshift range 
considered) in the cooling run. We attribute this 'positive' 
feedback effect to the metal enrichment active only in the 
former (we have already demonstrated that the gas accretion 
rates are very similar in both runs). Increased metallicity 
means more efficient cooling and this combined with the 
wind which can spread the metals, means a greater mass of 
gas will cool and condense enough to form stars. In fact, we 
find that the mass outflow rate is significantly lower than 
the rate of input of ejecta and entrained ISM into the SNe 
blastwaves, suggesting much of this material is recycled for 
next generation star formation. 

(iv) Star formation is globally very efficient in both the 
cooling and feedback runs. Even though in the star forma- 
tion implementation we set the efficiency to 0.01 when the 
gas density in a cell goes over 10 5 atoms/cc, we measure av- 
erage efficiencies of ~ 0.2 within 0.1r v ir- This indicates that 
star formation at high redshift proceeds in dense clumps that 
represent a larger fraction of the total gas mass of the galaxy 
than they do at low redshift. The global efficiency is well 
correlated with the average inflow rate of clumpy gas (i.e. 
gas that is condensed in satellites/subhaloes) (see Fig. 1131 
bottom) i.e. with the occurrence of mergers or strong inter- 
actions between neighbouring galaxies. When mergers occur 
the ISM is compressed and/or turbulence triggers increased 
fragmentation, allowing a larger fraction of the gas to reach 
high enough densities for star formation, resulting in higher 
global star formation efficiencies. 



8 APPENDIX: GAS PHYSICS 

8.1 Cooling and heating 

The cooling and heating of gas is included in RAMSES 
by adding an additional source term to the energy equa- 
tion. Above T = 10 6 K the cooling is dominated by 
bremsstrahlung radiation whereas in the range 10 4 < T < 
10 6 K, collisional and ionization excitation and recombina- 
tion processes dominate and can cool the gas to T — 10 4 
K. To cool below T — 10 4 K, molecular cooling and/or fine 
structure metal cooling need to be implemented; in the sim- 
ulations presented here, cooling down to IK is possible via 
metal line cooling. We note that collisional ionization equi- 
librium is adopted. 

Altho ugh it does not m a tter f or the current paper, we 
turn on a lHaardt fe Madaul (|l996h like UV radiation field 
(which varies with z, but is uniform in space) instanta- 
neously at z = 8.5 to model the reionization of the Universe. 

8.2 Tracking the metals 

The metallicity of the gas is followed as a passive scalar 
and is advected with the flow of gas. For both simulations 
with and without SNe feedback, we force the pristine gas to 
have a constant floor metallicity Z — 0.001Z©. This value 
provides a fairly good approximation to the real H2 cooling 
rates at work when the Universe is still devoid of heavy 
elements (Yohan Dubois, priv. comm.). For the simulation 
with SNe feedback, this floor is soon superseded by higher 
values because of the injection of heavy elements produced in 
the explosions of massive stars. When a SN explosion occurs, 



a fixed percentage of newly synthesized metals in the star 
particle (the 'yield', which we set to 0.1) is returned to the 
gas when creating the Sedov blastwave (see the section on 
supernova feedback). 



8.3 Star formation 

T he implementation of star form ation in RAMSES is outlined 
in iRasera fe Tevssierl (|200d ) and iDubois fe Tevssierl (|2008T ) 
and we recall the main details here. Given observational evi- 
dence for a universal relation between the gas surface density 
and the star for mation rate over many orders of magnitude 
at low redshift (|Kennicuttl fl998). we assume star formation 
proceeds according to a Schmidt law (a volumetric represen- 
tation of the Kennicutt law) . Gas is therefore converted into 
stars in cells in which the gas density, p, exceeds a threshold 
value, po, at a rate given by, 




where t* is the star formation scale, given by, 

/ \ -1/2 

,.=*,(£) (7) 

and is proportional to the free-fall time. 

To set the star formation parameters, we refer to mod- 
ern observations which have shown that star formation oc- 
curs at a constant efficiency of ~ 1 — 2% of gas per free-fall 
time over 5 orders of m agnitude in density, up to 10 5 cm -3 
(jKrumholz fe T an 2007). While it is not possible to achieve 
stellar densities in our simulations, we set the density thresh- 
old as high as our spatial resolution allows us, i.e. pj, the 
density at which the Jeans length associated with gas at 
100 K becomes equal to Ax m i n ~ 0.5 pc. This allows us 
to identify potential sites of star formation as realistically 
as possible. In the simulations described in this paper this 
corresponds to a density threshold at the top of the range 
currently probed by observations, po = 10 5 cm -3 , and a com- 
patible star formation efficiency of 0.01. Note that such a 
high density threshold for star formation allows us to by- 
pass the issue of having to decide which fraction of the gas 
is molecular and thus eligible for star formation as virtu- 
ally all the gas is in molecular form at such densities. Since 
the free-fall time for gas at po = 10 5 cm -3 is ~ 0.2 Myr, we 
then set to = 20Myr to obtain our chosen efficiency. Finally, 
in order to completely avoid artificial gas fragmentation on 
scales of the order of the Jeans length the following poly- 
tropic equation of state is also enforced: 

T=To (iT~ 1 (8) 

where To is set to 100K to prevent the Jeans length from 
becoming shorter than Ax m i n and 70 to 4/3 to stabilize the 
gas against gravity. 

When the conditions for star formation are met in a 
cell, the number of star particles to be created, N, is drawn 
from a Poisson distribution, 
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P(N) 



where 



pAx 3 

TTl* mil 



At 

t, 



(9) 



(10) 



where Ax 3 is the cell volume and At is the timestep eval- 
uated at the current refinement level. The minimum star 
mass, ra* )m in, is related to the resolution of the simulation 

by, 



1 + VSN + ??W 



(11) 



where Ax min is the volume of a cell on the finest refinement 
level of the simulation, and ?7sn and rjw are dimensionless 
parameters that account for the mass loss of stars through 
SNe explosions and the entrained gas mass respectively. De- 
tails of ?7sn and r]w are given in section 18.41 below. In our 
simulations, ra*, m i n ~ 167M , but note that this is reduced 
to ra*, m in ~ 77M© in the feedback run, since ?7sn and r]w 
are non-zero. 

In order to limit the number of star particles spawned 
in a simulation, the N particles which the Poisson law tells 
us should be created in a given cell are then fused together 
into one, so that each star particle mass is ra* = iVra*, m i n . 
To preserve numerical stability, no more than 90 % of the 
gas in any given cell is turned int o stars in one times t ep (th is 
is the same threshold as used bv lDubois fe Tevssierl §008)). 



8.4 Supernova feedback 

The scheme fo r SNe feedback emp l oyed in our feedback run 
is presented in iDubois &; Tevssierl (|2008h . When a star par- 
ticle is formed the amount of gas depletion can be broken 
down into the mass locked up in stars permanently, ra*, and 
the mass which will be ejected by the supernova, m e j, where 



= ra*(?7sN + rjw) 



(12) 



?7sn is the fraction of the total mass of stars formed which is 
recycled and becomes supernova ejecta. We set ?7sn — 0.106, 
i.e. the standard value for a Salpeter IMF truncated between 
0.1 and 100 M©. ?7w is a free parameter, the mass loading 
factor which quantifies the mass of gas surrounding the star 
that will be swept away by the shock wave which propagates 
outwards after the explosion, on scales that we cannot re- 
solve; this additional mass is added to the star particle when 
it forms. Setting r]w = 1, as we do, ensures that angular mo- 
mentum is conserved when we enforce t he Sedov blast wave 
solutio n described below for the gas (see lDubois &; Tevssierl 
(2008), Appendix A for details). Note that this means that 
the same amount of gas is entrained in the explosion as that 
whic h is lo c ked in to stars. Observations of galactic outflows 
bv iMartinl (|l999h show that the removal of gas from the 
galactic disc occurs at a rate that is comparable to the rate 
of star formation, which somewhat supports our choice of 
parameter, even though these measurements are performed 



on a much larger scale than that at which we enforce our 
Sedov solution ( ~1 pc). 

Star particles undergo SNe explosions lOMyr after their 
formation, corresponding to the typical lifetime of a 10M© 
Type II supernova progenitor. The total energy released, 
-dejecta, is set using the typical energy, Esn = 10 51 ergs, and 
mass, Msn = 10M©, of a Type II supernova, 



dejecta 



m» 

nsN M^ EsN 



(13) 



Half of this energy is kinetic and produces a Sedov blast 
wave with an initial radius of 2 cells on the highest level of 
refinement (~ 1 pc). The other half of the energy is thermal 
and heats the surrounding gas. Note that at these very high 
resolutions, the results are not very sensitive to the exact 
form (kinetic or thermal) in which the energy is deposited 
(Y. Dubois priv. comm). 
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